! ***************************************
/clear
/filn, E54  
/title, FEA of connecting zone of nozzle to cylinder
/units,si                        !ùʵλ
! ********* 趨*********
Rci=1000                         ! Ͳڰ뾶
tc=30                            ! Ͳ
Rco=Rci+tc                       ! Ͳ뾶
Lc=4000                          ! Ͳ峤
Rno=530                          ! ӹ뾶
tn=15                            ! ӹܺ
Rni=Rno-tn                       ! ӹڰ뾶
Li=193                           ! ӹ쳤
Ln=500                           ! ӹ쳤
rr1=30                           ! Բǰ뾶
rr2=15                           ! ڲԲǰ뾶
pi=6                             ! ѹ
pc=pi*Rci**2/(Rco**2-Rci**2)     ! Ͳ˲ƽغ
!****************ǰ***************************
/prep7
et,1,95                          ! 嵥Ԫ    
et,2,45                          ! 嵥Ԫ         
mp,ex,1,2e5                      ! ϵĵģ
mp,nuxy,1,0.3                    ! ϵĲɱ
TB,BISO,1,1,2,  
TBTEMP,0
TBDATA,,235,,,,,
!****************ģ***************************
cylind,Rco,Rci,0,-Lc/2,90,270,   ! Ͳ
wpoff,0,0,-Lc/2                  ! -ZƶLc/2
wprot,0,90,                      ! yzת90
cylind,Rno,Rni,-Ln-Rci-tc,-Rci+Li,90,180,  ! ɽӹ
vovlap,all                       ! overlap
vsel,s,,,7                       ! ѡͲ 
*afun,deg                        ! 趨ǶȺеλΪǶ
ang1=2*nint(asin(Rno/Rci))       ! ӹиǶ
wprot,0,0,-90+ang1               ! תϵ
vsbw,all                         ! иͲ
afillt,21,12,rr1                 ! ͲӹԲ  
afillt,23,35,rr2                 ! ͲڱӹԲ
afillt,14,25,rr2                 ! ¸Բ
afillt,13,19,rr1                 ! ϸԲ
alls
askin,91,64                      ! ݽӹԲڽӹϵĽ
vsba,4,13                        ! иӹ
askin,83,72                      ! ݽӹڹԲڽӹϵĽ
vsba,5,18                        ! иӹ
askin,67,75                      ! ԲǱ߽Ƥи24
vsba,2,24                        ! иͲ
wprot,0,,90-ang1                 ! תϵ
wprot,0,90,0                     ! תϵ
wpoff,0,0,-2*Rni                 ! ƶϵӹиλ
vsbw,all                         ! иͲ
VDELE,       3, , ,1             ! ɾͲϿʧȥĲ
ADELE,      25, , ,1             ! ɾ 
ADELE,      35, , ,1             ! ɾ
al,33,65,5                       ! ԲĶ
al,49,73,40                      ! ԲĶ
al,66,7,37                       ! ԲĶ
al,41,74,50                      ! ԲĶ
asel,s,,,31                      ! ѡԲ
asel,a,,,38,39
asel,a,,,4,6,2
va,all                           ! Բ
asel,s,,,21                      ! ѡԲ
asel,a,,,41,42
asel,a,,,5,9,4
va,all                           ! Բ
vsel,s,,,1,5,2                   ! ѡ
vsel,a,,,8,9
vsel,a,,,6
vadd,all                         ! ϲ
alls
aadd,4,43,29,26,5,51             ! ϲ
aadd,50,9,52,28,45,6             ! ϲ
LESIZE,67, , ,80, , , , ,1       ! ָԲܷʷ
LESIZE,5, , ,5, , , , ,1         ! ָԲʷ
vsweep,10                        ! ɨʷָ
alls
type,2
MSHAPE,0,3D                      ! 趨άӳʷ
MSHKEY,1 
VMESH,7                          ! ӳʷֽӹ϶
VMESH,4                          ! ӳʷֽӹ϶
accat,56,65                      ! Ͳ忪
LESIZE,108, , ,40,0.5, , , ,1    ! ָͲӹʷ
VMESH,13                         ! ɨʷָ
vsweep,11                        ! ɨʷͲ
vsweep,12
vsweep,2
alls
nummrg,all                        ! ϲͬ
numcmp,all                        ! ѹ
fini                             
!*******************************************
/solu
asel,s,loc,z,0                    ! ѡͲ
SFA,all,1,PRES,-pc                ! ʩӶƽغ
asel,s,loc,x,0                    ! ѡԳ
asel,a,loc,z,-Lc/2                ! ѡԳ
DA,all,SYMM                       ! ʩӶԳԼ
asel,s,loc,y,Ln+Rci+Tc            ! ѡӹܶ
DA,all,UY,                        ! Լλ
asel,s,,,38                       ! ѡڱ
asel,a,,,34
asel,a,,,9
asel,a,,,41
asel,a,,,12
asel,a,,,10
asel,a,,,7
asel,a,,,22
asel,a,,,25
asel,a,,,15
asel,a,,,27
asel,a,,,23
cm,ainner,area
outres,all,all
pred,on
KBC,0
autots,on
NLGEOM, on 
nsubst,1000
time,1                                 ! غɲ1
SFA,all,1,PRES,pi                      ! ʩѹ
alls
solve
fini
!*******************************************
/post1
PLNSOL,S,INT,0,1                  ! ʾӦͼ
fini